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Multi-valent (MV) battery architectures based on pairing a Mg metal anode with a high-voltage 
3 V) intercalation cathode offer a realistic design pathway toward significantly surpassing the energy 
storage performance of traditional Li-ion based batteries, but there are currently only few electrolyte 
systems that support reversible Mg deposition. Using both static first-principles calculations and 
ab initio molecular dynamics, we perform a comprehensive adsorption study of several salt and 
solvent species at the interface of Mg metal with an electrolyte of Mg^+ and CU dissolved in liquid 
tetrahydrofuran (THF). Our findings not only provide a picture of the stable species at the interface, 
but also explain how this system can support reversible Mg deposition and as such we provide insights 
in how to design other electrolytes for Mg plating and stripping. The active depositing species are 
identified to be (MgCl)^ monomers coordinated by THF, which exhibit preferential adsorption on 
Mg compared to possible passivating species (such as THF solvent or neutral MgCU complexes). 

Upon deposition, the energy to desolvate these adsorbed complexes and facilitate charge-transfer is 
shown to be small (~ 61 - 46.2 kJ moU^ to remove 3 THF from the strongest adsorbing complex), 
and the stable orientations of the adsorbed but desolvated (MgCl)^ complexes appear favorable 
for charge-transfer. Finally, observations of Mg-Cl dissociation at the Mg surface at very low THF 
coordinations (0 and 1) suggest that deleterious Cl incorporation in the anode may occur upon 
plating. In the stripping process, this is beneficial by further facilitating the Mg removal reaction. 


I. INTRODUCTION 

Since their commercial introduction nearly twenty-five 
years ago, rechargeable Li-ion batteries have performed 
admirably as the linchpin technology enabling today’s 
mobile electronics industry, currently powering hundreds 
of millions of laptops, cameras, and phones worldwide.® 
However, after years of continued improvement in perfor¬ 
mance we are approaching a fundamental limit of what 
can be accomplished with the current rocking chair Li- 
ion battery technology platform.® Further improvements 
in energy density (the quantity of energy that can be 
stored per charge or discharge either by mass or volume) 
without increasing the overall cell cost are still required, 
not only to lengthen the per-charge battery life of mobile 
devices, but also to facilitate the widespread adoption of 
electrified vehicles and grid-scale renewable energy stor¬ 
age. In short, a disruptive innovation in electrochemical 
energy storage can bring about tremendous technological 
and societal change— it has the realistic potential to dis¬ 
place significant CO 2 emissions from two of the largest 
global industrial contributors in transport and electricity 
generation. 

One of the most promising approaches is to build on 
the traditional robust three-component Li-ion battery 
architecture (intercalation or metal anode, non-aqueous 
electrolyte, and intercalation cathode), but with a multi¬ 
valent chemistry such as where multiple elec¬ 

trons are transported per ion (i.e. 2 electrons for each 
Mg^+ compared to 1 for Li+). This feature increases the 
energy density per intercalated ion in an electrode, but 


equally important, along with the switch to a Mg-based 
chemistry comes the potential to use a Mg metal anode 
as opposed to an intercalation structure (^ 700 Ah/L for 
Li in graphite compared to 3830 Ah/L for metallic Mg). 
For Li the use of metallic anodes especially challenging 
due to safety reasons.® 

The path to fully functioning and cost-competitive Mg- 
ion batteries begins with addressing a few but complex 
scientific questions, and perhaps the most pressing is to 
develop a robust understanding of the atomistic mecha¬ 
nism of reversible plating (or deposition) and stripping 
(or dissolution) of Mg at the anode/electrolyte interface 
during battery operation. To date, reversible Mg plating 
with low over-potential and reasonable anodic stability 
has been achieved in practice with only a specific class 
of electrolytes, namely organic or inorganic magnesium 
aluminum chloride salts (magnesium-chloro complexes) 
dissolved in ethereal solvents^H^. Aurbach and collab¬ 
orators developed the first operational Mg ion full cells 
using electrolytes based on Grignard reagents (e.g. butyl- 
magnesium chloride) ,® and improved upon their perfor¬ 
mance by switching to an in-situ prepared magnesium 
organo-haloaluminate, which can achieve high coulombic 
efficiency during stripp ing and deposition and anodic sta¬ 
bility up to 3.0 Recently, Muldoon and co-workers 

further improved upon air-sensitivity and anodic stabil¬ 
ity by designing a non-nucleophilic electrolyte containing 
a Hauser base, such as hexamethyldisilazide magnesium 
chloride (HMDSMgGl) and AlCls.^^ Another functioning 
electrolyte, though displaying limiting anodic stability, 
can be achieved using a combination of Mg(BH )4 and 
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LiBH 4 in diglyme as demonstrated by Shao et De¬ 
spite this substantial progress, some aspects of the elec¬ 
trolytes are in need of improvement, particularly the an¬ 
odic stability, safety and compatibility with other cell 
componentsJHlln contrast, most other non-aqueous elec¬ 
trolytes, including the Mg analogues to the common Li- 
ion battery electrolytes, lead to immediate electrolyte de¬ 
composition on the Mg anode, which passivates the metal 
surface preventing further electrochemical react ions P 

Much effort has been devoted toward characterizing 
the precise sequence of atomic-scale processes that make 
up the overall Mg deposition/dissolution process in ethe¬ 
real solvents. First, extensive experimental and compu¬ 
tational approaches have been used to characterize the 
equilibrium solvation structure of Mg^+ in the bulk elec¬ 
trolyte, which defines the thermodynamic start (end) 
point of deposition (dissolution), and several possible co¬ 
ordinating complexes involving com binations of Mg^+, 
Cl“, and THF have been identified.These 
findings have subsequently been used to inform phe¬ 
nomenological models that describe the general sequence 
of ion desolvation, adsorption on the anode surface, 
charge-transfer, and metal incorporation that make up 
the overall deposition process (the sequence in reverse 
describes the dissolution process) see Ref. |TT] and refer¬ 
ences therein. Only if all of these individual processes 
are kinetically unobstructed an Mg anode can operate 
reversibly and efficiently, e.g. with low overpotential and 
high coloumbic efficiency. 

However, the complete understanding at the atomic 
scale of the processes of electrochemical plating and strip¬ 
ping at the Mg anode/electrolyte interface, which would 
enable rational design of the next-generation Mg-based 
electrolytes, is still lacking. In recent work. Doe et 
showed that purely inorganic-based electrolytes solutions 
(i.e. without organometallic moieties), also referred as the 
“MACC” electrolyte, can exhibit reversible Mg electro¬ 
deposition with anodic stability up to 3.1 V vs. Mg. 
Hence, the MACC electrolyte system can serve as an ex¬ 
cellent prototype system to understand Mg deposition 
and stripping in the presence of Mg-Cl complexes. Here, 
we present an in-depth first-principles study aimed at 
characterizing the relevant chemical structures and com¬ 
plexes that form at the interface of a Mg metal anode in 
contact with an electrolyte solution comprised of Mg^+ 
and CD dissolved in an ether-based solvent (tetrahydro- 
furan or THF). By performing static first-principles ad¬ 
sorption calculations combined with ab initio molecular 
dynamics (AIMD) simulations, we are able to separately 
probe the interaction of the THF solvent in contact with 
the anode surface, as well as the full electrolyte/anode 
interaction by incorporating both Mg^+ and CD. Com¬ 
bined, these calculations reveal important features about 
magnesium-chloro complexes and their ability to func¬ 
tion in practice: first, cyclic ether-based solvents are con¬ 
firmed to be chemically inert against the Mg anode, while 
polymer-like ethers display comparatively more passivat¬ 
ing behavior on the Mg anode; and second, our calcula¬ 


tions reveal the existence of several (MgCl)+ coordina¬ 
tion complexes similar in energy but solvated by different 
numbers of THF molecules that define a kinetically facile 
Mg^+ desolvation process during deposition. 

The computational and theoretical strategy used to 
gain insights into the species at the Mg anode electrode 
and Mg desolvation/dissolution is non-specific and can be 
readily applied to a variety of important problems such 
as electroforming of metals, corrosion of metals and al¬ 
loys, medicinal chemistry, and catalysis. We emphasize 
that to obtain physically meaningful results, we consider 
explicitly the interaction of non-aqueous liquids with a 
solid phase, which is commonly approximated by a va- 
por/solid interface in other studies. 

II. METHODOLOGY 

Absorption calculations and ab initio molecular dy¬ 
namics (AIMD) simulations were performed withi n the 
DFT approximation, usi ng th e vdW-DF functionaP^^^ 
as implemented in VASP.^^^^ The total energy was sam¬ 
pled on a well-converged 3x3x1 A;-point grid together 
with projector augmented-wave theor}®! and a 520 eV 
plane-wave cutoff, and forces on atoms were converged 
within 1x10“^ eV A“^. For the different adsorption ge¬ 
ometries of Figure and Figure the adsorption energy 
(AE) is defined as: 

AF/ = F/gurf.+Ads. (s) F/gurf. (s) F/Ads. (1) : (1) 

where the Dsurf.+Ads.(s) is the total energy of the species 
formed by the Mg-surface and the adsorbing molecule 
(solvent or salt) in the solid state, Egurf. (s) the total en¬ 
ergy of the surface in the solid state, and EAds.O) 
total energy of the n adsorbing species. The DAds.(l) is 
conventionally approximated by the total energy of its 
gas reference, but to address the specific problem of the 
electrode in contact with a liquid electrolyte it is more 
practical to refer to the liquid state. To compute the liq¬ 
uid reference DAds. (1) we proceed by isolating low-energy 
snapshots from AIMD runs of both salt (MgCD or MgCl 2 
in THF solvent) and solvent (THF) and optimize their 
geometries. The average total energy provides a good es¬ 
timate of the liquid reference and incorporates both sta¬ 
bilizing 1st solvation shell effects and vdW-effects from 
the solvent. 

As the coordination of (MgCl)+ by THF in bulk is 
still debated, two different liquid (MgCl)+ bulk references 
were computed: i) four coordinated Mg, (MgCl)+-3THF 
as determined by Wan et al^ and ii) six-coordinated 
Mg, (MgCl)+-5THF as proposed by previous experimen¬ 
tal work.isnraini 

For the AIMD calculations the cutoff was reduced to 
400 eV and the total energy was only sampled at the 
F-point. Using a sampling of 1 fs, a production period 
of 15 ps was preceded by an equilibration period of 20 
ps within the canonical ensembl e (NV T) at 300 K, us¬ 
ing the Nose-Hoover thermostatThe slabs for the 
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adsorption calculations were cleaved from the fully re¬ 
laxed geometry of the Mg bulk (with a = 3.202 A and 
c = 5.139 A) and MgO (with a = 4.201 A). Mg (0001) 
and MgO (100) surfaces were modeled using a 5 layer 
periodic slab. Specifically for the adsorption of THF 
and poly-THF a4x4(a = b = 12.808 A) slab with 
80 Mg atoms was employed, whereas a2x2 (a = b = 
8.931 A) slab with 90 atoms for MgO. For the adsorp¬ 
tion of (MgCl)+ and MgCl 2 complexes, a larger surface 
with size 5x5 (a = 6 = 16.010 A) and 125 Mg atoms 
is needed. A well-converged vacuum of 35 A was in¬ 
serted between surface images, and the adsorbing species 
were always adsorbed on both surfaces to minimize artifi¬ 
cial electric dipoles through the non-periodic direction of 
the slab. The atomic positions of the two most exposed 
layers of the slabs and the adsorbates were optimized. 
All the geometries employed in the AIMD simulations 
were created by inserting a number of THF molecules at 
the experimental density (0.889 g/cm^) in the vacuum 
space between slabs to simulate liquid, resulting in unit 
cells with ^ 550 atoms. To accelerate computationally 
costly AIMD simulations, we performed an intermediate 
1 ns pre-equilibration NVT classical molecular dynam¬ 
ics simulation using the C HARMM potential as imple¬ 
mented in while keeping the coordinates 

of the surface and adsorbate atoms fixed as obtained from 
preliminary adsorption calculations. The concentration 
of the MgCD salt in the AIMD is 0.30 M and falls in 
the range of the experimental concentrations 0.25 - 0.40 
M explored by Aurbach and collaborators .^^Diffusivities 
were extracted from the mean squared displacements of 
the molecules at the surface from the AIMD simulations. 




FIG. 1. (Color online) (a) Stable adsorption sites (TOP, 
HLW, FLAT, and BRG) for THF on Mg (0001), and corre¬ 
sponding (b) adsorption energies AE (in kJ mol“^ referenced 
to the THF liquid state) with snapshots of their fully relaxed 
THF configurations (inset), (c) Variation of the adsorption 
energy with THF surface coverage O(THF)/Mg. 


III. RESULTS 
A. THF on Mg 

Adsorption energies computed from first-principles cal¬ 
culations are used to quantify the chemical interaction of 
a solvent tetrahydrofuran (THF) molecule with the Mg 
metal anode. From the standpoint of electrostatic inter¬ 
action, the THF molecule is expected to preferentially 
align with the ethereal oxygen oriented toward the ex¬ 
posed Mg sites at the surface. Previous experimentaP^ 
and theoreticaP^ work has demonstrated that Mg depo¬ 
sition and growth occur preferentially on the Mg (0001) 
facet at low current densities (^ 0.5 — 1.0 mA cm“^), 
hence the analysis is restricted to the Mg (0001) surface. 
Figure and[^ show four plausible adsorption sites at 
the Mg (0001) surface for a THF molecule labeled TOP, 
BRG, HLW, and FLAT, and their adsorption energies 
(AE). In this manuscript we refer strictly to the hep hol¬ 
low site HLW as indicated by the triangle in Fig. as 
opposed to the less stable fee hollow site, which shows to 
have consistently weaker adsorption energies. 

The calculated adsorption energies of Figure (ref¬ 
erenced to liquid THF) fall into a narrow range be¬ 


tween -10 and -17 kJ mol“^ with negligible variation 
between distinct surface sites and are consistent with 
the inert nature of ethereal molecul es in the pres¬ 
ence of a highly electropositive metal.l^im^ Negative 
AEs indicate that spontaneous adsorption is preferred. 
Referencing to the THF gas state would make the ad¬ 
sorption energies more negative (see Supplementary In¬ 
formation Figure SI) by adding the gas to liquid bond¬ 
ing energy (^ -46 kJ mol“^). The lower AE (^ -17.2 
kJ mol“^) of the FLAT configuration originates from the 
THF molecule’s orientation at the surface, favoring not 
only the 0-Mg(surface) interaction (of the top config¬ 
uration) but also additional weaker C-Mg(surface) and 
H-Mg(surface) bonds. THF molecules initially adsorbed 
in the BRIDGE sites are found to relax to the FLAT 
configuration (see Figure [^). 

A more realistic picture of THF adsorption, consistent 
with experimental conditions, is achieved by considering 
the effect of THF coverage. Figure shows the varia¬ 
tion of the adsorption energies for increasing number of 
adsorbed THF. The interaction of THF with the surface 
noticeably decreases with increasing THF coverage, and 
the adsorption energies shift to further positive values 
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FIG. 2. (Color online) PDF, g(r) of ethereal oxygen of 
THF and exposed Mg (0001) surface atoms (black lines), 
and surface coverage (dashed red line) for (a) THF, (b) 1,4- 
dihydroxy butane, and (c) DME (1,2-dimethoxyethane) 
from AIMD simulations. Insets are representative snapshots 
of THF and poly-THF fragments at the Mg surface. Red 
and blue shades highlight the lower bound of PDFs and the 
2'^^ solvation shell, respectively. 


consistent with the occurrence of THF-THF lateral elec¬ 
trostatic repulsions. As a comparison, we also report AEs 
that do not include van der Waals interactions (marked 
PBE in see Figure H), which provide a helpful upper- 
bound limit for the adsorption energies. The entropic 
contributions originating from the THF adsorption from 
the liquid will further increase the adsorption free ener¬ 
gies to more positive values, displacing the THF molecule 
from the Mg-anode and ensuring that ether molecules do 
not passivate the surface. 

To test the findings provided by the adsorption cal¬ 
culations and further resolve the dynamic environment 
of the anode/solvent interface, AIMD simulations within 
the canonical ensemble were performed at room temper¬ 
ature (300 K) and at the experimental density of THF 
(0.889 g/cm^).l^ Figure (black line) shows the pair 
distribution function, PDF or g(r), of the exposed Mg 
surface sites versus the THF oxygen atoms, Mgsurf.-O. 

The two dominating peaks observed at ^ 2 and 2.5 A 
show the co-existence of two THF configurations at the 
Mg(OOOl) surface, FLAT (peak at ^ 2.0 A) and TOP 
(peak at ^ 2.5 A), which also have the most negative 
adsorption energies, as seen in Figure [^. By integrat¬ 
ing the PDF profiles of Figure (red line) we conclude 
that a maximum of four THF molecules (2 in FLAT 
and 2 in TOP conformation) per 16 Mg exposed sur¬ 
face sites can be simultaneously adsorbed (dashed red 
line in Figure |^), consistent with the adsorption results 
of Figure [^. The diffusivity of THF molecules at the 
surface extracted from their mean squared displacement 


observed in the AIMD simulation is (6.1 ±1.5) x 10“^ 
cm^ s“^. Physically, the mean squared displacement of 
a molecule at the metal surface is related to the strength 
of adsorption— if the interaction is very strong (attrac¬ 
tive), molecules near the surface remains predominantly 
in the adsorbed state rather than in the liquid, resulting 
in lower surface diffusivity. 


B. poly-THF on Mg 

It is well understood that either electrophilic attack by 
protons on the ethereal carbon (self-polymerization) or 
catalytic interaction with an electropositive metal (such 
as A1 in AICI 3 or Mg) can initiate the THF ri ng openi ng, 
followed by polymerization forming poly-THF ! ^ Q We 
therefore also investigate the possible adsorption of THF 
products of polymerization such as poly-tetrahydrofuran 
(poly-THF) at the anode/electrolyte interface shown in 
the inset of Figure Here, two approximations are 
used to describe poly-THF: i) by cleaving THF at the 
ethereal bonds and saturating by hydrogen atoms, form¬ 
ing the 1,4-dihydroxybutane (see Figure]^), and ii) by 
1,2-dimethoxyethane also known as glyme (DME, see 
Figure [^). While the acidic nature of the terminal 
hydrogens in 1,4-dihydroxybutane may enhance the ad¬ 
sorption energy with the Mg-surface compared to linear 
ethers (e.g. glyme, tetraglyme, nglyme), it nonethe¬ 
less provides an upper bound to the adsorption energies 
of a variety of THF polymerization products that might 
be present in the cell. The adsorption energies of 1,4- 
dihydroxybutane (referencing to a dilute solution of 0.4 
M 1,4-dihydroxybutane-THF in liquid THF) for differ¬ 
ent configurations were always found to be more neg¬ 
ative (-32.0 ± 5 kJ mol“^) than those of cyclic THF. 
Examination of the Mgsurf.-O PDF reveals that the oxy¬ 
gen atoms of 1,4-dihydroxybutane are always nearer to 
the Mg-surface (broad peak from 1.7 to 2.7 A) in com¬ 
parison to cyclic THF molecules (see Figure |^, b), 
likely due to more geometric degrees of freedom ac¬ 
cessible by 1,4-dihydroxybutane that can allow for im¬ 
proved binding for a linear molecule on the anode sur¬ 
face. Measurements of the diffusion coefficients (from 
AIMD simulations) further indicate a smaller mobility 
of 1,4-dihydroxybutane at the surface, ^ (9.5 ± 1.6) 
xl 0 “^ cm^ s“^, as compared to weakly physisorbed and 
more mobile THF molecules, ^ (6.1 ± 1.5) x 10“^ cm^ 
s“^. However, at the Mg(OOOl) surface DME is found 
to bind less strongly (-10.5 ± 1 kJ mol“^) than 1,4- 
dihydroxybutane (-32.0 ± 5 kJ mol“^), and comparable 
to THF (-17.2 kJ mol“^), confirming the inert nature of 
this molecule. Although the similarity of the Mgsurf.-O 
PDFs of DME and THF (see Fig. and Fig. [^) at¬ 
tests comparable binding nature for these molecules, the 
minimum distance between the ethereal oxygen of DME 
and Mg atoms at the surface falls at slightly larger dis¬ 
tances ^e peak at 2.1 A) than in THF (see peak at 2.0 A 
in Fig. ^), which can be explained by the larger steric 
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hindrance of DME. 

Recently, Barile et al^ have observed that the addi¬ 
tion of poly-THF in the MACC electrolyte increases the 
deposition overpotentials to 500 mV indicating that poly- 
THF may passivate more the Mg-surface. We conclude 
that hydroxyl-terminated polymer molecules (e.g. 1,4- 

dihydroxybutane) are more prone to passivate the Mg- 
anode surface (than linear-ethers and THF), thereby neg¬ 
atively affecting the battery performances during strip¬ 
ping and deposition. Furthermore, the chelating nature 
of linear ether with long chains (e.g. glyme, diglyme, 
tetraglyme, and poly-TH F) is e xpected to inhibit Mg^+ 
ion delivery at the surface. 


C. THF on MgO 

The irreversible reaction of Mg with ubiquitous water 
forms resistant passivation layers that suppress all Mg 
electrochemical activity and is one of the limiting factors 
in the development of Mg-ion batteries. Investigations 
of Mg corrosion have clarified the bilayer structure of 
the passivating film, with the inner pa rt comp rised of 
MgO and an external layer of Mg(OH) 2 p^^^^^^In order 
to clarify the effect of surface passivation on solvent ad¬ 
sorption, we calculated the adsorption energies for THF 
on the MgO (100) termination,!^ which is found to grow 
epitaxially on the Mg(0001) facet. Analogous to previ¬ 
ous sections we use static DFT adsorption calculations 
to address the interaction of THF with MgO(lOO), and 
verify these findings with auxiliary ab initio MD simu¬ 
lations. We anticipate that the ethereal oxygen of THF 
will interact preferentially with Mg atoms at the surface 
(as also seen by THF on Mg metal) while the THF’s H 
atoms will maximize H-bonding with surface oxygens. 

Figure 1^ and[^ show four plausible adsorption sites 
at the MgO (100) surface for a THF molecule labeled 
TOP, FLAT, BRG and BRG-H and their respective ad¬ 
sorption energies (AE). BRG-H represents a THF ad¬ 
sorption arrangement that maximize H-bonding between 
THF hydrogen atoms and MgO oxygen species at the 
interface (as seen in Fig. [^). TOP, FLAT, and BRG 
configurations are the same as described for adsorption 
on Mg. 

Surprisingly, the adsorption energies reported in 
Fig. (referenced to liquid THF) are found remarkably 
close to AEs computed for THF on Mg(OOOl), varying 
between -17 and -18 kJ mol“^. By close examination 
of each THF adsorption configuration and relative AEs, 
we observe that the BRG configuration relaxes to TOP, 
and similarly the FLAT configuration relaxes to BRG-H, 
which maximizes the interactions with the surface oxygen 
atoms through hydrogen bonds (see snapshots of Fig. [^). 
On the basis of these results we cannot observe a sub¬ 
stantial difference between the adsorption properties of 
THF on MgO (100) and Mg (0001) surfaces and suggest 
that further theoretical and experimental investigations 
are needed. 



FIG. 3. (Color online) (a) Stable adsorption sites (TOP, 
FLAT, BRG and BRG-H) for THF on MgO(lOO), and corre¬ 
sponding (b) adsorption energies AE (in kJ moP^ referenced 
to the THF liquid state) with snapshots of their fully relaxed 
THF configurations (inset). THF ethereal oxygen in blue for 
clarity. 


AIMD simulations are performed to investigate the dy¬ 
namic properties of THF molecules at the MgO inter¬ 
face with similar settings described in previous sections. 
Figure]^ and[^ (black line) show the pair distribution 
functions, PDF or g(r), of the exposed Mg and O surface 
sites versus the THF oxygen atoms, Mgsurf.-OxHF and 
Osurf.-OxHF, respectively. 

Similar to Fig. the Mgsurf.-OxHF PDF of Fig. 
shows two major peaks at ^ 2 and 2.5 A, emphasizing 
the co-existence of two THF adsorption configurations 
at the MgO(lOO) surface: BRG-H (peak at ^ 2.0 A) and 
TOP (peak at ^ 2.5 A). The remarkable similarity of the 
Mgsurf.-O PDFs on Mg(OOOl) and MgO(lOO) (seen com¬ 
paring Fig. and Fig.[^) suggests that PDF analysis is 
not sufficient to distinguish whether THF molecules are 
adsorbed on active or passivated Mg anodes. The PDF 
Osurf.-OxHF shown in Fig.[^ is more informative and can 
discriminate between adsorption on active or passivated 
surfaces, with a single peak between 2.7 and 4 A charac¬ 
teristic to absorption on MgO. Mgsurf.-OxHF PDF also 
provides information about the THF coverage on MgO 
(see red line in Fig. |^), concluding that a maximum of 
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FIG. 4. (Color online) (a) PDF, g(r) of ethereal oxygen of 
THF and exposed Mg surface atoms and (b) O exposed sur¬ 
face atoms on MgO(lOO) (black lines), respectively. Coverage 
(dashed red line). Insets are representative snapshots of THF 
fragments at the MgO surface. Red and blue shades high¬ 
light the lower bound of PDFs and the 2'^^ solvation shell, 
respectively 

two THF molecules (1 in FLAT and 1 in BRG-H con¬ 
formation) can be simultaneously adsorbed per 12 Mg 
exposed surface sites (dashed red line in Fig. [^). 

D. Magnesium-chloro Salt and THF solvent 

By first establishing that ethereal solvent molecules 
do not significantly interact with the Mg anode sur¬ 
face, we turn our attention to characterizing the an¬ 
ode/electrolyte interface by including the magnesium 
chloride salt species (solvated by THF). The atomic- 
scale interactions between salt and solvent that de¬ 
fine the relevant reacting complexes and species in 
the bulk electrolyte can be quite complex as demon¬ 
strated in seminal work by Aurbach and collabora¬ 
tors, who characterized intricate equilibria between 
magnesium-chloro clusters and A1C1^“ in the form 
of “simpler” monomers [Mg(/i-Cl)-5THF] + [(AlCl 4 )]“, 
dimers [Mg 2 (/ i-Cl) 3 - 6 TH F] + [(AlCl 4 )]~, and larger poly¬ 
meric speciesA recent theoretical work by 
Wan et alW clarified the first solvation shell of the 
monomer and dimer magnesium organo-chloro species 
in the bulk electrolyte. Their AIMD calculations sug¬ 
gest that the (MgCl)+ monomer is always coordinated 
by three THFs. However, it is still debated whether the 
(MgCl)+ monomer is actively present in the bulk elec¬ 
trolyte or only formed in proximity of the electrode as 
the product of the dimer decomposition.^^ The struc¬ 
ture of the organo-chloroaluminate salts (C 2 H 5 MgCl- 
[(C 2 H 5 ) 2 A 1 C 1]2 in THF, DCC) in proximity of the Mg- 
anode was investigated experimentally using XANES 


spectroscopy by Benzmayza et revealing that only 
the cation monomer [Mg(/i-Cl)'5THF]+ can approach the 
surface (as opposed to the more complex and bulkier 
dimer and trimer), thus representing th e active species 
during the Mg-deposition. Previous worlPESHIIIll on the 
organo and inorgano chloroaluminate electrolytes have 
elucidated the complex thermodynamic equilibria be¬ 
tween monomer, dimer and larger multimeric species, 
which might ultimately affect the concentrations of the 
(MgCl)+ near the anode electrode. Even though more 
complex Mg-carriers (e.g. dimer and larger multimeric 
units) could also be present at the surface, we set to 
investigate only the monomer (MgCl)+ as the simplest 
possible vehicle to transport Mg^+ ions from the elec¬ 
trolyte in proximity of the anode.^ In this work, we use 
these findings to inform the construction and preparation 
of our first-principles simulations, focusing on the inter¬ 
action of (MgCl)^ monomers coordinated by THE near 
the anode surface. 

Using a similar approach to the one used to study the 
anode/solvent interface, we first probe the adsorption en¬ 
ergies of possible magnesium-chloro complexes, and then 
perform AIMD simulations to explore more precisely the 
dynamic nature of the anode/electrolyte salt interface. 
Envisioning the overall Mg deposition process, for exam¬ 
ple, beginning with a THE-solvated (MgCl)+ complex 
in the bulk electrolyte and ending as metallic Mg incor¬ 
porated in the anode, the deposited Mg must in some 
sequence shed its sheath of THE molecules. In Eigurej^ 
the interfacial structure and corresponding adsorption 
energies of several charged magnesium-chloro (MgCl) + 
complexes (solvated by 5 through 0 THE molecules and 
in different orientations) adsorbed on a Mg(OOOl) sur¬ 
face are shown. Throughout the manuscript we adopt an 
abbreviated nomenclature to discuss magnesium-chloro 
(MgCl)+ complexes adsorbed at the Mg surface, namely 
the first number indicates the number of THE ligands 
followed by the orientation (or site) adopted by the ad¬ 
sorbing (MgCl)+ salt. Eor example, 4-SIDE means that 
(MgCl)+ is coordinated by 4-THE molecules and lays 
horizontally on the Mg surface as shown in the snapshots 
of Eigurej^ Two different liquid references were used to 
compute the AEs of adsorbed magnesium-chloro salts at 
the Mg surface: i) bulk (MgCl)+-3THE reference,^ and 
ii) with a bulk (MgCl)+- 5 THF,En 3 indicated respec¬ 
tively by the upper and lower bounds of the bars in Fig¬ 
ure The difference in absorption energies from these 
two reference states is small, and can be attributed to 
the slightly more stable (MgCl)+-3THF reference in the 
reference calculation. 

Upon first inspection, the magnitude of the energies 
for salt adsorption (in Figure compared to THF ad¬ 
sorptions (values shown in Figure and text above) 
confirms that the THF solvent weakly interacts with the 
Mg metal surface while magnesium-chloro salt complexes 
readily adsorb, in good agreement with experimental ob¬ 
servations by Benzmayza et a/.^ an d the reversible plat¬ 
ing observed in this system.^l ^ * ^ * ^^ * ^^ *In all cases. Cl aligns 
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FIG. 5. (Color online) stable configurations of (MgCl)^ solvated by varying No. of THF (from 5 to 0) adsorbed at the Mg 
(0001) surface, and corresponding adsorption energies (AE in kJ moP^). Upper and lower bounds indicated by the bars on 
AEs are from using 3-THE and 5-THE bulk liquid reference states, respectively. Black circles and color-coded arrows indicate 
THE removal. 


toward the Mg surface, with adsorption on the HLW site 
preferred over the TOP site as shown in Figure (also 
see Figure S2 in Supplementary Information), because in 
comparison to TOP, HLW and SIDE oriented adsorption 
(shown in Figure]^ maximizes the interaction of the Cl 
ion with the Mg surface sites as well as that of some 
dangling C and H atoms from the THF molecules with 
the Mg surface. Of all the charged magnesium-chloro 
complexes, (MgCl)+ coordinated by 5 THF in the hol¬ 
low site exhibits the strongest adsorption energy (-84.7 
± 10.5 kJ moP^), and in general the strength of the ad¬ 
sorption reduces with sequential THF removal, as seen 
in the trend of adsorption energies of Figure Also, as 
THF is successively removed from the salt complex, the 
most stable orientation relative to the surface shifts from 
HLW to SIDE. Both trends can be understood by con¬ 


sidering the role of the Mg surface in stabilizing known 
preferred Mg^+ coordination geometries.ll^ For 5-HLW, 
Mg^+ can remain stably in octahedral coordination as 
seen in Figure but as more THF are removed; Mg^+ 
is increasingly exposed to unstable coordinations. For 
example, 4-HLW has a more positive adsorption energy 
in comparison. Adopting the SIDE orientation allows 
Mg^+ to be additionally coordinated by the Mg surface, 
which can stabilize the salt complex when there are too 
few THFs to “ideally” solvate Mg^+ as seen in Figure 
where the SIDE configuration becomes the most stable 
for 3 THFs or fewer. As even more THF molecules are 
removed, the adsorption of the salt complexes is contin¬ 
ually weakened, but a notable exception is from 3-SIDE 
to 2-SIDE (as seen in Figure where there is an in¬ 
crease in adsorption strength likely due to less steric hin- 








drance between THF molecules which can be relieved by 
removing the most weakly bound THF. The adsorption 
of (MgCl)+-lTHF,1-SIDE, and bare (MgCl)+, 0-SIDE 
on the Mg(0001) surface were found to be very unfavor¬ 
able (+96.8 and +31.1 kJ mol“^ respectively), indicating 
that these species are likely too unstable to play a signif¬ 
icant role in the reversible Mg deposition process. 

As mentioned earlier, the (MgCl)+ monomer sur¬ 
rounded by a solvent coordination shell of 4 or 5 THE 
molecules adsorbs strongly to the Mg surface (as seen 
in Eigurej^, but these complexes also preserve the in¬ 
tegrity of the weak Mg-Cl bond. In addition to electron 
transfer, completion of the Mg deposition process also re¬ 
quires the eventual dissociation of Mg from Cl“ at the Mg 
surface, which necessitates breaking the solvation shell. 
The differences in adsorption energies shown in Eigurej^ 
can also be considered as the energy to (de) solvate the 
(MgCl)+ complex by sequential (removal) addition of 
THE molecules, and initial desolvation appears facile 
with only a cumulative ^ 46.2 - 61.0 kJ mol“^ required to 
remove 3 THEs from the stably adsorbed (MgCl)+-5THE 
complex. To investigate this observation further under 
more realistic conditions comparable to experiment, we 
perform room temperature AIMD simulations with dif¬ 
ferent (MgCl)+ complexes adsorbed on the Mg surface in 
a liquid THE reservoir. Eigurej^ shows the pair distribu¬ 
tion function between the Mg atom of the (Mg-Cl) + salt 
and the oxygen atoms of the THE ligands (black lines) 
and the respective coordination numbers (red lines) for 
three different salt adsorption configurations, 5-HLW, 4- 
SIDE and 2-SIDE. In Eigurej^the average Mg-O(THE) 
distance is centered around 1.9 - 2.2 A, falling in the 
typical Mg-0 bond length range (2.1 - 2.2 A). The coor¬ 
dination number measured at the center of the g(r) peaks 
(see red dashed lines in Eigurej^as a guide for the eye) is 
4, 3, and 2 for 5-HLW, 4-SIDE and 2-SIDE, respectively. 

Surprisingly, in some cases the number of THEs in the 
first coordination shell of the (MgCl)+ salt spontaneously 
decreases by one THE unit after few picoseconds of equi¬ 
libration time in the AIMD simulation. Eor example, 
when beginning the AIMD simulation with either 5-HLW 
or 4-HLW absorbed on the Mg surface, we observe as the 
final result 4-HLW and 3-SIDE, respectively, as seen in 
the dashed red lines in Eigurej^ and Eigurej^. Likely, 
this is due to additional THE-THE interactions between 
the ligands of the adsorbed complex and bulk liquid in the 
AIMD simulation (unaccounted for in the static adsorp¬ 
tion calculations which are performed in vacuum) and 
the effect of entropy, which stabilizes THE in the liquid 
compared to the adsorbed state. This evidence further 
supports Mg ion desolvation as a facile process. As seen 
in Eigure the trend of spontaneously shedding THE 
stops in the 2-SIDE configuration, as the significant in¬ 
crease in desolvation energy required to further remove 
THE (as seen in Eigure cannot be offset by the finite 
temperature and additional solvent interactions. 

The large positive energy increases when removing ad¬ 
ditional THE i.e. from 2-SIDE to 1-SIDE to 0-SIDE, 


suggests that, the electron transfer likely occurs before, 
when there are still some surrounding THE molecules 
(likely 2), which further underscores the necessity of con¬ 
sidering the stabilizing coordination effect of the sol¬ 
vent ion with an explicit coordination shell of THE. 
Erom close inspection of the relaxed geometries of low 
and zero-coordinated THE complexes (such as adsorbed 
1-SIDE and 0-SIDE, see Eigure [^, we observe incipi¬ 
ent dissociation of the Mg-Cl bond (this also occurs for 
some of the MgCl 2 configurations (see Eigure S5 in the 
Supplementary Information), which further suggests that 
magnesium-chloro complexes with few THE molecules 
are the relevant species involved in charge transfer. Eor- 
tunately, the SIDE orientation is most energetically fa¬ 
vorable at lower THE coordinations (as shown in Eig¬ 
ure |^, which allows for the depositing Mg^+ to become 
better exposed to the anode surface, which should better 
facilitate ensuing charge transfer. 

Einally, we also consider the adsorption on the Mg- 
surface of the charge-neutral MgCl 2 species that may 
precipitate at the anode interface, for example as undis¬ 
solved/unreacted MgCl 2 during the preparation and pos¬ 
sibly conditioning^ of the magnesium-chloro electrolyte. 
Previous AIMD simulations by Wan et a/.^ have shown 
that MgCl 2 in bulk THE solvent remains coordinated by 
2 THE molecules (MgCl2-2THE), which we adopt as the 
liquid reference to compute the adsorption energies. In 
general, the AEs are always found to be very positive for 
all the models considered in this study i.e. bare MgCl 2 
(174.8 kJ mol-i), MgCl2-2THE (25.0 kJ moD^), and 
MgCl2-3THE (28.9 kJ mol“^), indicating that MgCl 2 ad¬ 
sorption is not competitive with the (MgCl)+ salts at 
equilibrium (see Eigure S6 in Supplementary Informa¬ 
tion). 


IV. DISCUSSION 

To gain practical insights into the design of electrolytes 
capable of supporting reversible Mg deposition in Mg-ion 
batteries, we performed a comprehensive first-principles 
calculations based adsorption study of the relevant salt 
and solvent species at the interface of a Mg metal sur¬ 
face in contact with a solution of Mg^+ and CD dissolved 
in THE solvent using static adsorption calculations and 
ab initio molecular dynamics simulations. Although our 
model relies on a few but important assumptions, namely 
a Mg metal anode free of defects, impurities, our results 
reveal important features of the system consistent with 
reversible Mg deposition and other experimental obser¬ 
vations. 

Erom static adsorption calculations of the an¬ 
ode/solvent system, THE exhibits weak interaction with 
Mg metal, which should leave the surface unpassivated 
and therefore accessible to salt adsorption. Eurthermore, 
the adsorption calculations of (MgCl)+ complexes co¬ 
ordinated by several THE molecules indicate that sol¬ 
vated salt complexes rather than solvent species are 
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models (a) 5-HLW, (b) 4-SIDE and (c) 2-SIDE (see Figure |^. Insets show the structural evolution of the salts at the surface. 
Red dashed lines as guide for the eye. Black circles indicate THF removal. Mg-Cl PDFs in Figure S4 of the Supplementary 
Information. 


strongly preferred at the surface, which is confirmed un¬ 
ambiguously in the AIMD simulations. Linear poly-THF 
molecules, however, reported to be the product of AICI 3 
catalyzed THF polymerization,^^! might show larger in¬ 
teraction with the anode surface in comparison to THF 
and linear ethers. 

In this investigation we have clarified that oxygen im¬ 
purities (in form of MgO) do not affect the adsorption 
properties of th e solvent, but as demonstrated by previ¬ 
ous studie^^^^ traces of other impurities such as Cl and 
A1 can potentially react with salt species and the anode 
participating to the passivation of the electrode. 

Considering the Mg deposition process, once a 
(MgCl)+ complex solvated by some (4 or 5) THF is ad¬ 
sorbed at the Mg surface, breaking the THF solvation 
shell to facilitate electron transfer is shown to be initially 
facile, requiring only ^61-46 kJ mol“^ to remove 3 
THF ligands from 5-HLW, for example. Simultaneously, 
the adsorbed complex adopts a conformation that brings 
the adsorbed Mg^+ ion nearer to the Mg surface, rotat¬ 
ing from HLW to SIDE orientation, which should benefit 
electron transfer and eventual Mg incorporation (by max¬ 
imizing the orbitals overlap between Mg^+ ions and Mg^ 
atoms of the electrode). The energy to remove additional 
solvent molecules abruptly increases once there are only 
two remaining THF molecules, and (MgCl)^ complexes 


solvated by 1 and 0 THF molecules do not even bind to 
the Mg surface, which strongly suggests that (MgCl) + 
coordinated by 2 THF molecules is the relevant complex 
involved in charge transfer, an important finding requir¬ 
ing further experimental verification. 

The electrochemical activity of Cl as an additive to 
facilitate the electron transfer during multi-valent metal 
deposition is well known for Cu deposition, demonstrated 
from both a theoretical and experimental standpoint by 
Nagy et alW^. The presence of Cl“ accelerates the first 
electron transfer from Cu^+ to Cu+ in the multistep re¬ 
duction of copper, and the overlap of Cl“ and Cu^+ or¬ 
bitals eases the electron transfer mechanism by isolat¬ 
ing the plating Cu^+ ions from competitive solvent in¬ 
teraction. In that model, it is crucial that the Cl-Cu^+ 
complex pre-adsorb at the surface, and indications of the 
same behavior in this system (preferential adsorption of 
(MgCl)+ complexes) are overall in agreement with ex¬ 
perimental observations of reversible Mg deposition. 

Although we did not explicitly consider the effect of an 
applied electrochemical potential to the system, the ad¬ 
sorption calculations in this work also shed some light 
on deposition events that occur post charge transfer. 
It is well established that atomic and ionic Cl readily 
adsorb in the hollow sites of Mg metal surfaces with 
highly exothermic adsorption energies (-241.8 kJ mol“^ 
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calculated in this work and in agreement with previous 
wor ]Pi), and the observed stability of isolated Cl species 
in our calculations suggests that lingering Cl adsorption 
on Mg surface after charge transfer may affect deposi¬ 
tion kinetics. Indeed, recent SEM, energy-dispersive X- 
ray spectroscopy and XPS measurement s reve al traces of 
Cl (in addition to Al) at the Mg surface.^I^^ In order to 
keep depositing pure Mg metal at the anode, residual Cl 
must be continually removed which can contribute to the 
plating overpotential, and we speculate that the addition 
of AICI 3 species to the electrolyte solution may improve 
de posit ion kinetics by “scrubbing” the Mg anode free of 
Cipi^ Specifically, added AICI 3 can act as a shuttle for 
Cl“ from the anode to the cathode during Mg plating 
by shifting the AICI 3 + Cl“ ^ AICI 4 reaction equilib¬ 
rium, as AICI 3 reacts with accumulated Cl“ at the anode 
surface, the AIC 4 product is driven toward the cathode. 
Continual removal of chloride at the anode interface is a 
source of overpotential upon deposition, but in the strip¬ 
ping process, the additional availability of Cl“ should 
promote the process. During the dissolution process, the 
applied potential drives the charged (MgCl)+ complexes 
to the cathode, which consequently increases the local 
concentration of the neutral MgCl 2 complexes at the an¬ 
ode, and as we observe, these spontaneously dissociate 
into (MgCl)+ and Cl“, thus creating more carriers to 
drive the reaction. Overall, this creates a source of asym¬ 
metry between the deposition and stripping process. 

V. CONCLUSIONS 

In conclusion, we have set to resolve the atomistic 
mechanisms that allow for reversible Mg electrodeposi¬ 
tion in non-aqueous electrolytes by performing a com¬ 
prehensive first-principles adsorption study at the in¬ 
terface of a Mg metal anode in contact with an elec¬ 
trolyte solution consisting of Mg^+ and Cl“ dissolved 
in THE (tetrahydrofuran) solvent. Erom an analysis 


of both solvent and salt adsorption from static first- 
principles and ab initio molecular dynamics calculations, 
we have gained a number of insights as to why this narrow 
class of electrolytes (magnesium-chloro salts dissolved in 
ethereal solvents) can function in practice: first, neither 
cyclic ether solvent molecules nor neutral salt complexes 
(MgCl 2 based complexes) exhibit surface-passivating be¬ 
havior (although some polymeric THE shows more favor¬ 
able surface interaction in comparison to cyclic THE); 
second, the strongest adsorbing species’ at the Mg inter¬ 
face are charged (MgCl)+ complexes which are also the 
active species involved in charge-transfer; third, the en¬ 
ergy to (de)solvate (MgCl)+ is minimal; and finally, the 
stable orientations of de-solvated (MgCl)+ complexes are 
also favorable for charge transfer. In future work, we plan 
to directly investigate the effect of an applied potential 
on the interfacial properties considered in this study with 
the goal of directly characterizing the electron-transfer 
mechanism upon Mg deposition and stripping. 
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